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ABSTRACT 

The dynamical signatures of the interaction between galaxies in clusters and the 
intracluster medium (ICM) can potentially yield significant information about the 
structure and dynamical history of clusters. To develop our understanding of this phe- 
nomenon we present results from numerical modelling of the galaxy/ICM interaction, 
as the galaxy moves through the cluster. The simulations have been performed for a 
broad range, of ICM temperatures {kTd = 1, 4 and 8 keV), representative of poor 
clusters or groups through to rich clusters. 

There are several dynamical features that can be identified in these simulations; 
for supersonic galaxy motion, a leading bow-shock is present, and also a weak gravita- 
tionally focussed wake or tail behind the galaxy (analogous to Bondi-Hoyle accretion). 
For galaxies with higher mass-replenishment rates and a denser interstellar medium 
(ISM), the dominant feature is a dense ram-pressure stripped tail. In line with other 
simulations, we find that the ICM/galaxy ISM interaction can result in complex time- 
dependent dynamics, with ram-pressure stripping occurring in an episodic manner. 

In order to facilitate this comparison between the observational consequences of 
dynamical studies and X-ray observations we have calculated synthetic X-ray fiux 
and hardness maps from these simulations. These calculations predict that the ram- 
pressure stripped tail will usually be the most visible feature, though in nearby galaxies 
the bow-shock proceeding the galaxy should also be apparent in deeper X-ray observa- 
tions. We briefiy discuss these results and compare with X-ray observations of galaxies 
where there is evidence of such interactions. 

Key words: galaxies: clusters: kinematics and dynamics - galaxies: interactions - 
intergalactic medium - X-rays: galaxies - galaxies: clusters: Virgo 



1 INTRODUCTION 

Most galaxies are not isolated (TuUy 1987) but are mem- 
bers of groups or clusters which are gravitationally bound, 
though not necessarily virialised. Galaxies in such systems 
may experience a number of types of interaction with their 
environment: galaxy-galaxy encounters, which may trigger 
starbursts and can culminate in galaxy merging (Schweizer 
1982), interaction with the potential of the system as a 
whole, which may for example, result in tidal stripping of 
loosely bound material from galaxy outskirts (Mamon 1987) , 
and finally hydrodynamical interaction between the inter- 
stellar medium (ISM) in a galaxy, and any intergalactic gas 
trapped in the potential well of the galaxy system. 

It is the last of these processes which is the subject of the 
present study. Interaction with the surrounding intracluster 
medium (ICM) may have a significant impact on the gas 
within a galaxy whenever it sweeps through a substantial 



mass of hot gas in its motion within a group or cluster. 
In practice, the hot ICM dominates the baryon content of 
galaxy clusters (White & Fabian 1995) and is also a major 
baryonic component within many galaxy groups (Mulchaey 
et al. 1996, Ponman et al. 1996). 

The interaction between galaxies and the ICM should 
generate characteristic structures: at low galaxy velocities 
it will produce only a mild displacement of the outer re- 
gions of any gaseous galaxy halo, at higher velocities a bow 
shock should be generated ahead of the galaxy, and a wake 
of enhanced gas density may trail behind it. 

These processes are important for three main reasons: 
first, gas may be removed from and/or accreted by the 
galaxy in such interactions with its gaseous environment, 
with important implications for subsequent galaxy evolu- 
tion. Second, such processes will be responsible for some or 
all of the metals which are observed today within the ICM 
in clusters and groups, and may also affect the energetics of 
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the ICM (Ponman, Cannon & Navarro 1999) - as such they 
are important for an understanding of the structure and evo- 
lution of galaxy systems. Finally, structures such as galaxy 
wakes provide information about the direction of motion of 
galaxies on the plane of the sky, which has in general bccu 
unavailable previously (Merrifield 1998). Given a good un- 
derstanding of the astrophysics of the interaction, one may 
also be able to deduce the magnitude of the galaxy's total 
velocity (allowing a full dynamical solution) and also to con- 
strain the total mass of galaxies, including any dark halo. 

Given the velocities of galaxies in groups and clusters, 
the temperature of the gas in shocks and wakes will be 
~ 10^ K, so the observational evidence of such processes 
is best sought at X-ray energies. To date shock and wake 
structures have only been observed in a few instances (de- 
tailed below), due to the limited sensitivity of existing X- 
ray telescopes. However, the improvements available from 
the upcoming X-ray observatories, Chandra and XMM will 
enable detailed studies of structures around galaxies within 
clusters to be conducted. The aim of the present paper is to 
prepare for this, by exploring the Eistrophysics of the galaxy- 
ICM interaction using numerical hydrodynamic simulations, 
and in particular to investigate the observable properties of 
such interactions in the X-ray band. 

The astrophysics is quite complex - several different 
mechanisms may be at work in any particular galaxy, and 
these will potentially have different signatures. Examples of 
such mechanisms include: i) a bow-shock, if the galaxy is 
moving supersonically through the ICM then there is the 
likelihood of a bow-shock preceeding the galaxy; ii) a gravi- 
tationally focussed wake, the mass and motion of the galaxy 
will tend to concentrate material behind it, in an analo- 
gous manner to Bondi-Hoyle accretion in binary systems 
(c./. Hunt 1971); iii) ram pressure stripping of gas in the 
galaxy by the ICM will also result in a tail of enhanced den- 
sity behind the galaxy. This feature, in the case of M86, has 
also been referred to as an X-ray plume (Rangarajan et al. 
1995). 

Examples of these features can be seen later in simu- 
lations of the galaxy/ICM interaction. While a bow-shock 
will be a distinct and recognisable feature, distinguishing 
between wakes (formed by gravitational focussing) and tails 
formed by stripping might be difHcult on morphological 
grounds. In this paper we identify spectral as well as mor- 
phological signatures that can be used as a diagnostic of the 
physical mechanisms at work. 

Currently there is only a limited amount of evidence 
of bow-shocks, wakes and tails around galaxies, observed 
primarily in relatively local systems. Some examples are: 

(i) White et al. (1991) reported on EINSTEIN observa- 
tions of the Virgo galaxy M86, finding evidence for an X-ray 
plume, caused by ram pressure stripping as the galaxy moves 
through the cluster ICM. Rangarajan et al. (1995) analysed 
ROSAT observations of M86, finding that the X-ray plume, 
which they also interpret as being due to ram-pressure strip- 
ping, has a higher metallicity than either the cluster or M86 
itself. The plume temperature is kT = 0.8 keV, similar to 
the galaxy emission from M86, but substantially lower than 
the cluster. 

(ii) David et al. (1994) reported on ROSAT observations 
of a linear structure behind the dominant galaxy in the 



NGC 5044 group, they interpreted these results as being due 
to a "cooling wake" behind the galaxy, caused by motion of 
the galaxy with respect to the hot gas which follows the 
potential of the galaxy group. 

(iii) Irwin & Sarazin (1996) reported on ROSAT HRI ob- 
servations of the Virgo cluster elliptical galaxy NGC 4472, 
which shows a distorted X-ray morphology. Irwin & Sarazin 
(1996) suggested that this was strong evidence for a bow- 
shock like structure around the galaxy. 

(iv) Sakelliou, Merrifield & McHardy (1996) reported 
on ROSAT results on observations of the wide-angled tail 
galaxy 4C 34.16 which lies in a poor cluster. Radio maps 
show two tails bent into a wide C-shapo, and X-ray observa- 
tions show an enhancement between the tails. Sakelliou et al. 
(1996) interpret this as due to motion of the galaxy through 
the ICM, and that this motion is responsible for both the 
jet bending and the X-ray wake trailing the galaxy. 

(v) The Fornax elliptical galaxy NGC 1404 shows an X- 
ray tail pointing away from the cluster centre, indicating 
that galaxy is infalling (Jones et al. 1997). 

(vi) ROSAT observations of the Coma cluster found an 
elongated X-ray structure around the giant elliptical galaxy 
NGC 4839 (Dow & White 1995). NGC 4839 is probably the 
cD galaxy of a subcluster that is currently infalling into the 
main cluster. The X-ray tail of NGC 4839 is pointing away 
from the main cluster centre, supporting the infall idea. 

(vii) Hq observations of the edge-on spiral galaxy 
NGC 4388 by Veilleux et al. (1999) have found evidence of 
supersonic ram-pressure stripping. NGC 4388 is moving at 
~ 1500 km s~^ with respect to the ICM of the Virgo cluster, 
with a Mach number oi A4 — 3. The X-ray emission from 
NGC 4388 is extended in the plane of the galaxy (Matt et 
al. 1994) but there is no evidence of X-ray extensions asso- 
ciated with ram-pressure stripping (though as Veilleux et al. 
(1999) point out, the galaxy orientation is unfavourable). 

There have been a several papers considering the gas- 
dynamics of a galaxy moving through a cluster since the 
issue was considered by Gunn & Gott (1972). Shaviv & 
Salpeter (1982) reported on simulations, concentrating on 
the role of galaxy motion in ram-pressure stripping. Further 
simulations of the galaxy/ICM interaction have been pre- 
sented by Takeda, Nulsen & Fabian (1984), who presented 
results of the dynamics of a galaxy moving radially through 
a cluster, suffering periodic stripping in the higher density 
regions near the centre of a cluster. Other studies of ram- 
pressure stripping have been presented by Nulsen (1982) and 
Portnoy, Pistinner & Shaviv (1993), as well as the work of 
Abadi, Moore & Bower (1999), which concentrates on the ef- 
fect of ram-pressure stripping in clusters on the truncation of 
disks of spiral galaxies. The two most relevant papers to this 
work are those by Gaetz, Salpeter & Shaviv (1987) and Bal- 
sara, Livio & O'Dea (1994). Gaetz et al. (1987) presented a 
comprehensive study of ram-pressure stripping from galax- 
ies, and calculated several important parameters afi^ccting 
the efficiency of stripping. We shall use some of the param- 
eters defined by Gaetz et al. (1987) in this paper. Balsara 
et al. (1994) present further calculations, using more sophis- 
ticated rmmcrical techniques. Importantly, they found that 
the stripping process was much more complicated than seen 
in the calculations of Gaetz et al. (1987). When there was 
significant mass-replenishment of the ISM within the galaxy 
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the resulting dynamical structure was complex, with mate- 
rial being stripped from the galaxy in long streamers from 
the edge of the galaxy, with the stripping happening is a 
discontinuous, episodic manner. 

The calculations presented here represent an improve- 
ment over those of Gaetz et al. (1987) and Balsara et al. 
(1994) for the following reasons. The Gaetz et al. (1987) pa- 
per used a rather old numerical hydrodynamics code, which 
has substantial numerical viscosity. This has the effect of 
smearing out smaller scale features found in runs using more 
recent numerical hydrodynamic computer codes, such as the 
Piecewise Parabolic Method (PPM) . This is important as in 
many runs the stripping process occurs in discrete events, 
with material being stripped of the edges of the galaxy. It is 
likely that in the Gaetz et al. (1987) models these processes 
were not be resolved. Balsara et al. (1994) used a more so- 
phisticated code (PPM - essentially the same as used in this 
paper). However, Balsara et al. (1994) used a large value for 
the mass-replenishment rate (Mrep). A large value of Afrep 
has the advantage of highlighting the complex dynamics of 
the galaxy/ICM interaction, but care has to be taken not to 
make it unrealistically large (see Section 2.3.3). 

There are two further areas where this work is a signif- 
icant improvement. First, from the simulations we calculate 
synthetic flux and hardness maps, as well as surface bright- 
ness profiles. These synthetic results can then be compared 
directly to actual observations of structures around galax- 
ies in clusters. Such an approach will be necessary to fully 
understand the observed structures. Second, we perform a 
broad ranging parameter study, for parameters appropriate 
for a range of clusters, ranging from cool clusters or groups 
(with an ICM temperature of kT^i ~ 1 keV), up to hot clus- 
ters {kT^i ~ 8 keV). 

The main purposes of this paper are twofold. First, to 
identify observational signatures of the galaxy/ICM inter- 
action, which can be used in future detailed studies of indi- 
vidual objects, and second, to identify the circumstances in 
which galaxy/ICM interactions are most likely to be visible. 

The paper is organised as follows: in Section ^ we dis- 
cuss the hydrodynamical model used to investigate the phe- 
nomenon, including the model for the galaxy and the pa- 
rameters used in the simulations. In Section H results from 
the simulations are presented and compared with analytic 
models of the galaxy/ICM interaction, in Section ^ the ob- 
servable signatures of these simulations are discussed, and 
in Section H we discuss these results in the context of obser- 
vations and conclude. 



2 THE NUMERICAL MODEL 
2.1 Introduction 

To model the gas dynamics in and around a galaxy moving 
within a cluster we use the 2-D hydrodynamics code VH- 
1, which has been used extensively elsewhere in different 
contexts (Blondin et al. 1990; Stevens et al. 1992). The code 
uses the PPM algorithm and is good at resolving shocks 
in complex flows. It is basically very similar to that used by 
Balsara et al. (1994). The flow is assumed to be cylindrically 
symmetric. 

In these simulations the gravitational effect of the 



galaxy on the ICM are included, as well as replenishment 
of the galaxy ISM by stellar populations. The parameters 
for both the galaxy and the clusters they are embedded in 
are described below. 



2.2 The Model Galaxy Cluster 

One of the main goals of this paper is to investigate the ex- 
pected observational visibility of the galaxy/ICM interaction 
region. To do this we embed the results from the hydrody- 
namic simulations into a model cluster. We assume a spher- 
ical, isothermal cluster with a temperature Td and cluster 
sound speed a^i (calculated assuming a mean mass per par- 
ticle oi fi — 10"^* gm). Simulations for three different values 
of Tci ; a) Tci = 1 keV have been performed, representative of 
a poor cluster, such as Fornax (c./. 4C 34.16, Sakelliou et al. 
1996), b) Tci = 4 keV, representative of a more substantial 
cluster, such as Abell 1060 or AWM 7 (both of which have 
temperatures ~ 4 keV; Loewenstein & Mushotzky 1996), 
and c) Tci = 8 keV, indicative of a rich cluster such as Coma. 

We assume the following scaling relationships for the 
cluster X-ray luminosity [L^) and velocity dispersion {od, 
White, Jones & Forman 1997): 



L^erg s-i) =6.1 x lO''" (fcT,,)' 
a,,(kms'^) =371 (fcTrf)"-"' 



(1) 
(2) 



with kTci in keV. For the model clusters a /3 = 2/3 mass 
distribution has been assumed (see below). The X-ray lumi- 
nosity given by eqn. (^) is the total bolometric X-ray lumi- 
nosity, rather than the ROSAT luminosity. A cluster core 
radius of Rcore ~ 250 kpc and a distance of D = 50 Mpc is 
assumed for all clusters (c./. Jones & Forman 1984). 

Future studies of individual galaxies in clusters will be 
able to use values of Lx (and /3 and (Jd etc) appropriate for 
the cluster in question. However, these scaling relationships 
are adequate for a general study of the observational char- 
acteristics of structures around galaxies. For simplicity, the 
galaxy is assumed to be moving through material of con- 
stant density, at a distance of r = Rcore from the cluster 
centre. 

For P = 2/3, the surface brightness distribution S{w) 
of the model cluster as a function of projected distance from 
the cluster centre w is 



S{w) = 



So 



{1 + iw/Rcorerf/^ ' 



(3) 



where Rcore is the cluster core radius, and is the central 
surface brightness. This profile integrates to give the a total 
X-ray luminosity of Lx = Siv'^ SoRcore- Also, for P = 2/3, 
the gas density profile within the cluster pcl{r) varies as 

Pci(O) 



Pci{r) 



{1 + ir /Rcore)'') 



(4) 



with /5c; (0) the cluster central gas density. At r = Rcore the 
gas density pci{Rcore) = Pc;(0)/2. Details of the parameters 
for the model clusters are given in Table |^. 



2.3 The Model Galaxy 
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Table 1. The assumed parameters for the model clusters. The various symbols are defined in more depth in Section 



Parameter 


Symbol (Units) 


Cool Cluster 


Intermediate Cluster 


Hot Cluster 


Cluster X-ray luminosity 


Lx { erg s~^ ) 


6.1 X 10*2 


2.8 X 10** 


1.9 X 10*5 


Cluster core radius 


Rcore (kpc) 


250 


250 


250 


Velocity dispersion 


(Tci (km s~l ) 


371 


795 


1164 


Cluster temperature 


fcTrf (keV) 


1.0 


4.0 


8.0 


Cluster sound speed 


a^l (km s~^ ) 


520 


1040 


1460 


Cluster density at galaxy radius^ 


PcliRcore) (gm Cm-3 ) 


3.8 X 10-2* 


2.5 X 10-27 


5.8 X 10-27 


Central cluster surface brightness-^ 


So (erg cm~'^s~-'^arcmin~^ ) 


1.1 X 10-1* 


5.1 X 10-" 


3.4 X 10-12 



Notes for Table 1: 

1 In all simulations the galaxy is assumed to be at a radius of r = Rcore from the cluster centre. 

2 The quoted surface brightness is calculated assuming the cluster bolometric X-ray luminosity. The surface brightness in the ROSAT 
waveband will be correspondingly lower. 



2.3.1 Mass Distribution 

The parameters for the model galaxy used here are similar 
to those of Gaetz et al. (1987). The galaxy is assumed to be 
spherical, with the mass in the galaxy dominated by stars. 
The stellar density distribution within the galaxy is given 
by 



p*(r) = 



P.(0) 



(1 + W)2) 



(5) 



with Rc the core radius of the mass distribution {Rc = 
1.5 kpc in these simulations), and p*(0) is the central den- 
sity in the galaxy. We assume that the mass distribution of 
the galaxy extends out to a radius Rh ~ 32.2 kpc (Balsara 
et al. 1994). The total galaxy mass, Mgai, is then 



Mgai =4np,{0)Rl 



Rh 
Rc 



tan 



1 

\ Rc 



(6) 



In these simulations Mgai = 1.2 x 10^2 Mq , which then spec- 
ifies the galaxy central stellar density p*(0). The assumed 
parameters for the model galaxy are summarised in Table gI 



2.3.2 Galaxy Velocity 

For a galaxy cluster with a one-dimensional velocity disper- 
sion of ad the mean galaxy velocity will be Vmean = VSctci . 
The model galaxy is assumed to be moving at a velocity Vgai 
through the ICM, and for each of the three model clusters 
simulations with two values of Vgai, have been run, with 
Vgai = 1.0 or 1.5Vmcan. The relationship between a^i and 
cluster temperature (eqn.|2|) means that Vgai oc kT^(^^ . 

For these values of Vgai all the model galaxies are mov- 
ing supersonically, with Mach numbers M = 1.2 — 2.1. 
This range of values should be representative, though some 
galaxies of interest have even higher Mach numbers (c.f. 
NGC4388, Veilleux et al. 1999). The steeper dependence 
of Vgai on Tci compared to the sound speed, means that 
the galaxies in clusters with higher kTci have marginally 
higher Mach numbers. Details of the velocities for the dif- 
ferent models are in Table M. 



2.3.3 Mass replenishment within the galaxy 

In galaxies a combination of stellar mass loss and supernovae 
act to continuously replenish the gas within the galaxy. This 
mass-replenishment will affect the flow around the galaxy in 



several ways. First, it provides material that can be con- 
tinuously stripped and second, it acts as barrier (if dense 
enough) to the ICM, resulting in the formation of a bow- 
shock proceeding the galaxy. In these simulations the radial 
distribution of the mass replenishment is assumed to have 
basically the same form as that of the stellar distribution, 
that is 



(l + (r/ii,)2) 

for r < Rrcp, where Rrcp is the outer radius for mass- 
replenishment. The characteristic radius Ro is the mass re- 
plenishment half-radius - half of the mass-replenishment oc- 
curs inside a radius Rq. This prescription is equivalent to 
"Case 1" in Gaetz et al. (1987). 

This expression can be integrated to give an expression 
for the total mass-injection rate within the galaxy 



Mr 



■ 4-Rprep{0)Rl 



Rc 



tan 



1 / Rrcp \ 



(8) 



The value of prep(O) is determined by the assumed to- 
tal rate of mass injection into the galaxy Mrep. Assum- 
ing a galactic mass/light ratio of 20 (in solar units) gives 
L = 5 X 10^° 1/0 for the model galaxy. This mass/light ra- 
tio implies a modest dark matter halo. Faber & Gallagher 
(1976) estimated that the rate of mass-replenishment into an 
elliptical galaxy was 0.015 yr-^ x (10® Lq ), so the ex- 
pected mass-replenishment rate for a 1.2 x 10^2 Mq galaxy 
is ~ 1 Mq yr-i . As discussed by Gaetz et al. (1987), the 
relationship of Faber & Gallagher (1976) was estimated as- 
suming a single coeval populations of stars, with a current 
main sequence turn-off mass of 1 AIq . Younger populations 
of stars will enhance the mass-replenishment rate, and so the 
value of 1 Mq yT~^ should be considered as a lower limit. 
The results of Mathews (1989) suggest a value of Mrcp = 
5 Mq yr-^ as a reasonable upper limit, when there is mass- 
replenishment from a younger stellar population. Conse- 
quently, two sets of simulations with Mrcp = 1 Af© yr-^ 
and 5M0 yr-^ have been calculated. The temperature of 
the replenished gas is assumed to be Trcp ~ 1-87 x 10^ K, 
the same as in Balsara et al. (1994) and Gaetz et al. (1987). 

We note that Balsara et al. (1994) used a rather high 
value of the mass-replenishment rate (~ 15 Mq yr-^ ). Such 
a large value leads automatically to the formation of dense 
structures in the centre of the galaxy. However, even with 
lower values of Mrcp, under some circumstances, dense and 
complex structures are formed within the galaxy. 
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Table 2. Parameters for the model galaxy used in all the simulations. The symbols used are defined in more depth in Section 



Parameter 


Symbol 


Value 


Cralaxy mass 




l.Z X lU 


Galaxy central stellar density 




9.48 X 10^^^ g 


Galaxy core radius 


Rc 


1.5 kpc 


Galaxy outer radius 


Rh 


32.2 kpc 


Mass replenishment half-radius 


i?o 


8.7 kpc 


Mass replenishment outer radius 


Rrep 


15.5 kpc 


Mass replenishment temperature 


T 

^ rep 


1.87 X 10^ K 


Star formation density threshold 


Psf 


2.34 X 10-2* g 


Star formation temperature threshold 


Tsf 


2.0 X 10"' K 


Star formation timescale 




3.25 X 10'^ yr 



2.3.4 Star-formation 



2.5 Useful Parameters 



In several of the simulations presented in Section Q a consid- 
erable amount of cold, dense gas accumulates in the centre 
of the galaxy. We include the possibility of star-formation 
in the simulations, in an identical way to Balsara et al. 
(1994), which is essentially as a mass-sink in the centre of 
the galaxy. No feedback mechanisms (i.e. the formation of a 
young stellar population leading to an increase in the mass- 
replenishment rate) is included. Such refinements could be 
easily included in more sophisticated models. 

As implemented here the action of star-formation is to 
remove gas from the central regions of the galaxy. Mate- 
rial is assumed to form stars or molecular clouds and ef- 
fectively drops out of the hydrodynamic equations. When 
star-formation does occur the rate of change of gas density 
is given by; 



Psf 



4r. 



1 + tanh 



1 + tanh 



sf 



Q.lPsf 

- T 



O.IT, 



(9) 



In this equation psf and Tsf are the effective density and 
temperature thresholds above which and below which star- 
formation will occur. The star-formation time-scale r^f spec- 
ifies the time-scale on which the star-formation occurs. 
In these simulations psf — 2.4 x IQ-^* gm cm"^ , Tsf = 
20000K, and Tsf = 3 x lO'^ years. 

In our simulations star-formation only occurs in a sub- 
set of the runs (see Section and when it does occur it 
only happens in the innermost regions of the galaxy. Star- 
formation, as implemented here, does not have any substan- 
tial effect on the flow-properties outside of these inner re- 
gions, but it does prevent the densities in the central region 
of the galaxy growing, unchecked, to unphysical levels. 



2.4 Description of Simulations 

A total of 12 simulations have been run, four each for the 
cool, intermediate and hot clusters, with two values for the 
galaxy velocity and mass-replenishment rate. Details of the 
parameters for each model are given in Table ^. The simula- 
tions were all run on a 240 x 180 2-D cylindrical grid, with an 
extent of 2.4 x 10^^ cm by 1.8 x 10^^ cm (~ 80 x 60 kpc). At 
an assumed cluster distance of 50 Mpc each cell corresponds 
to an angular extent of 1.3 arcsec. 



Gaetz et al. (1987) introduced some parameters and 
timescales that are useful to characterise these simulations. 
The flow-crossing timescale is defined as 

2Ro 



'Tflow 



gal 



(10) 



For the range of velocities considered Tfiom is in the range 
(0.6 — 2.6) X 10^ yr. The characteristic orbital timescale 
within the galaxy is defined as 

2ttRo 



'^orb 



with 



( GMjRo) 
\ Ro 



1/2 



(11) 



(12) 



with M{Ro) the total mass-contained within a radius Ro. 
For the model galaxy M{Ro) = 2.6 x lO" , so that 
Vcircfi = 360 km s^^ , and Torb = 4.7 x 10^ yr. 

An important quantity in these simulations will be the 
fraction of mass retained by the galaxy in the gas-phase and 
the amount of gas ram-pressure stripped from the galaxy. 
Gaetz et al. (1987) parameterised the fraction of mass ram- 
pressure stripped from the galaxy using two parameters, u 
and uj, where 



TrRopclVgal 



M 

i\i rep 

with Pel the ICM density and 

Vgal 



LO = 4 



/ Vgal \ / Vesc,0 \ ' 

VlO^kms^V V500kms-V 



(13) 



(14) 



.103 km s"^ 

The escape velocity from the galaxy at a radius _Ro is 

W,o = (2|0(i?o)-<?i(oo)|)^''2 . (15) 

with the gravitational potential. For the parameters used 
here Vesc,o ~ 850 km s~^ (Table ^). 

The parameter is the ratio of the external to internal 
mass-flow rate into the galaxy, and u! a dimensionless mea- 
sure of galaxy velocity. With these two parameters Gaetz 
et al. (1987) defined a dimensionless stripping parameter 
^ = VLu^''^ that related directly to the fraction of mass- 
stripped from the galaxy (their Fig. 3). Consequently, the 
rate of ram-pressure stripping is oc PdVg^f. 

From the simulations we calculate the fraction of the 
replenished gas that is retained by the galaxy ($). We esti- 
mate the mass of gas Mgas within a radius Ro at two times 
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Density Plot — Model 2a Temperature Plot — Model 2a 




Distance (cm) Distance (cm) 



Figure 1. The density and temperature structure for two simulations that show the "Efficient Stripping" (ES) behaviour. Upper panels - 
results for Model 2a (4 keV cluster - low velocity and low mass-replenishment rate). Note the bow-shock structure, and also the relatively 
small mass density/temperature perturbations associated with the galaxy. Lower panels - results for Model 3d (8 keV cluster - high 
velocity and high mass-replenishment). The bow shock is visible and has a smaller opening angle due to the higher Mach number of the 
flow. 



ti and t2, with ti and t2 chosen to lie at times where the 
solution has settled down to a ste ady state (or at least a 
quasi-steady state - see Section 3.2). We then define $ as: 



$ : 



[Mgas{t2) - Mgas{tl)] + [M,/(t2) - Msf{tl)] 



(16) 



{t2 -tl)Mrep 

where Msf{t) is the total mass removed from the system by 
star-formation up to a time t. From their models, Gaetz et 
al. (1987) empirically found that $ could be represented by 
the following expression 



1 + 



6/2 



(17) 



with ^1/2 = 0.6 is the value of ^ at which $ = 0.5. Our 
results for $ will be discussed in Section ^ 

Other parameters from the simulations that will be 
useful for comparison with observations are the emission 
weighted temperatures of specified regions within the simu- 
lations. For all of the simulations we have calculated the 
emission weighted temperature for two regions. The first 
covers the galaxy itself, and comprises a sphere, radius i?o, 
centred on the galaxy (Tgai). The second covers the tail re- 
gion, and is a cylinder, radius -Ro, extending from a distance 
Ro behind the galaxy to the edge of the grid [Ttao)- The 
emission weighted temperature T is defined as 



r : 



(18) 



with n the gas number density and A(r) the emissivity of 
the gas. The summation is carried out over all cells within 
the region of interest. The values of T^i and Ttaii for each 
of the simulations are shown in Table pi although it should 
be noted that because some of the simulations never settle 
down to a steady state, the values of Ttaii and Tgai can vary 
with time. 



3 RESULTS 

In this section we discuss the basic dynamical features found 
in these simulations. In Section ^ we go on to discuss the 
expected observable features that may be seen at X-ray en- 
ergies. 

There are two basic dynamical modes of behaviour in 
these simulations (illustrated in Figs. |l| and The first is 
when the galaxy is completely stripped by the ICM (and 
the mass-retention fraction $ ~ 0), and there is no substan- 
tial build up of material at the galaxy centre. We term this 
flow behaviour the "Efficient Stripping" or ES mode. The 
second is when there is a substantial build up of material 
in the centre of the galaxy, and the stripping is less efficient 
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Figure 2. The density and temperature structure for two simulations for the "Mass-Retention" (MR) mode. Upper panels - results for 
Model lb (1 keV cluster, low velocity, high mass-replenishment rate). Note the strong density enhancements in the form of streamers in 
the tail region. There is a considerable amount of cool and dense gas at the centre of the galaxy. Lower panels - results for Model 2c 
(4 keV cluster, low velocity, high mass-replenishment rate). The bow shock is again clearly visible as well as a dense tail behind the 
galaxy caused by ram-pressure stripping. The tail is not as over dense as for the case in the top-panels, and is smaller in extent. 



($ ~ 0.1 — 0.6). We term this flow behaviour the "Mass- 
Retention" or MR mode. In the first case the dynamics are 
relatively simple, while for the second case the dynamics are 
complex, resuhing in a chaotic downstream region. Results 
from simulations showing both forms of behaviour are shown 
m Figs. |l| and |[ For each of the model runs listed in Table H 
the main dynamical mode is identified. As discussed later, 
however, the distinction between the modes is not always 
completely clear cut. 

3.1 The "EfHcient Stripping" Mode 

We illustrate the case of the "Efficient Stripping" mode by 
looking, in detail, at two representative simulations, one each 
for the intermediate and hot cluster models (Models 2a and 
3d). The density and temperature structure for both models 
are shown in Fig. |l| 

Consider Model 2a first. This is a lower Mach number 
simulation, with M = 1.3. The density and temperature 
structure, at a time 5 x 10* yr into the run, are shown in the 
top panels of Fig. |l]. The basic flow structure is relatively 
simple and the simulation quickly settles down quickly to a 
steady state. The main features that can be seen in Fig. ^ 
are as follows: 

(i) a prominent bow-shock in front of the galaxy, visible in 



both the density and temperature plot. There is an increase 
in density and temperature across the shock. The peak post- 
shock temperature is 5.7 x 10^ K, compared to the ambient 
ICM temperature of 4.6 x 10^ K (Fig. ^). This increase in 
temperature is that given by the Rankine-Hugoniot condi- 
tions for a shock transition with this Mach number. 

(ii) The gravitational acceleration of the ICM material by 
the galaxy results in a lower density region preceeding the 
galaxy. 

(iii) There is little material accumulated at the centre of 
the galaxy, with a peak density of p = 4 x lO"'^^ gm cm~^ 
(Fig. The density never gets high enough for star- 
formation to occur. In contrast, the peak density can be 
much higher in simulations with a higher Mrep- 

(iv) The galaxy is completely ram-pressure stripped, with 
the mass-retention fraction "1> ~ 0. 

(v) On account of the relatively small perturbations in 
density and temperature the values of Tgai and Ttaii are 
similar to the ambient cluster temperature. 

(vi) The ICM velocity fleld is not greatly affected by 
the galaxy. Material in front of the galaxy is gravitation- 
ally accelerated, then decelerated both by the shock and 
by the mass-loading of replenished material, which is swept 
straight out of the galaxy. There are no regions where there 
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Table 3. The assumed and derived parameters for the simulations. 



(km s-1 ) 



Ad rep 

{ Mq yr- 



M 



(K) 



^tail 

(K) 



Mode 





(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


Cool cluster models {kTd = 


1 keV = 1.2 X 


lO^K) 


















Model la 


640 


1.0 


1.2 


0.59 


0.33 


1.2 


X 


10^ 


1.1 X 


107 


MR 


Model lb 


960 


1.0 


1.9 


1.6 


0.31 


3.5 


X 


lO"* 


1.1 X 


107 


MR 


Model Ic 


640 


5.0 


1.2 


0.12 


0.89 


2.4 






1.1 X 


107 


MR 


Model Id 


960 


5.0 


1.9 


0.31 


0.75 


1.3 


X 


10* 


1.1 X 


107 


MR 


Intermediate 


cluster models (kT^i = 4 keV 


= 4.6 


X lO'^K) 
















Model 2a 


1380 


1.0 


1.3 


24.7 


0.0 


5.1 


X 


10^ 


5.1 X 


107 


ES 


Model 2b 


2070 


1.0 


2.0 


65.5 


0.0 


4.8 


X 


107 


5.2 X 


107 


ES 


Model 2c 


1380 


5.0 


1.3 


4.9 


0.03 


5.2 


X 


107 


2.9 X 


107 


MR 


Model 2d 


2070 


5.0 


2.0 


13.1 


0.0 


5.6 


X 


107 


4.9 X 


107 


ES 


Hot cluster models {kT^i = 


i keV = 9.3 x lO^K) 


















Model 3a 


2010 


1.0 


1.4 


141.5 


0.0 


9.5 


X 


107 


9.7 X 


107 


ES 


Model 3b 


3015 


1.0 


2.1 


374.5 


0.0 


9.4 


X 


107 


9.8 X 


107 


ES 


Model 3c 


2010 


5.0 


1.4 


28.3 


0.0 


9.9 


X 


107 


9.3 X 


107 


ES 


Model 3d 


3015 


5.0 


2.1 


74.9 


0.0 


9.8 


X 


107 


9.8 X 


107 


ES 



Notes on Table 
Column (1) 
Column (2) 
Column (3) 
Column (4) 
Column (5) 



Galaxy velocity with respect to the cluste r ICM (Section | 
Galaxy mass-replenishment rate (Section |j.3.3| ). 

The Mach number of the galaxy with respect t o th e undisturbed ICM. 
The ram-pressure stripping parameter (Section 2.5). 
The galaxy mass-retention fraction (eqn. 16h. 



Columns (6) and (7): The mean temperatures in the tail region and galaxy centre region (Section 2.5). 

Column (8): The mode of stripping for this model: ES - "Efficient Stripping"; MR - "Mass-Retention" (see Section ^ for an 
explanation) . 



is a region of downstream reverse flow back in the galaxy 
(c./.Fig. |). 

The results for Model 3d (shown in the lower panel of 
Fig. 1^), are essentially the same as for Model 2a. The re- 
sults are shown at a time 2.5 x 10* yr after the start of the 
simulation. In spite of the fact that the mass-replenishment 
rate for this simulation is higher than for Model 2a the ram- 
pressure of the cluster is still sufficient to completely strip 
the galaxy. From Fig. ^, for Model 3d the opening angle of 
the bow-shock is substantially smaller. This is due to the 
higher Mach number of the flow [M = 2.1) for Model 3d 
compared to = 1.3 for Model 2a. 

The expected observational signatures for these simula- 
tions will be discussed in Section 0. 



3.2 The "Mass-Retention" Mode 

In contrast to the relatively simple dynamics for the "Ef- 
ficient Stripping" cases described above, when the mass- 
replenishment is sufficient or the ram-pressure low enough, 
complex and time-dependent flow patterns result in the 
vicinity of the galaxy. In these cases, material accumulates in 
the central region of the galaxy and dense blobs and streams 
of material are stripped from the galaxy. In this situation we 
do expect to see observable consequences of the galaxy /ICM 
interaction. 

The MR mode is illustrated in Fig. ^ with results from 
two simulations (Model lb, top panels and Model 2c, bottom 



panels), and we will concentrate on Model lb initially. The 
density and temperature profiles shown in Fig. ^ are at a 
time 6.8 x 10* yr after the start of the simulation. The MR 
simulations typically do not settle down to a steady-state, as 
there is a radial "pumping" mode instability in the flow (as 
discovered and discussed by Balsara et al. 1994). However, 
the MR mode simulations do reach a broadly quasi-steady 
state configuration. 

As in the case of "Efficient Stripping" a large bow-shock 
region leading and enveloping the galaxy is apparent, caused 
by the interaction of ICM with the galaxy. As noted earlier 
we also a leading bow-shock in the case when efficient strip- 
ping occurs. In this case the opening angle of the bow-shock 
is broadened by the accumulated replenished material at 
the centre of the galaxy. Behind the bow-shock is a region 
of shock heated gas (though the temperature enhancements 
is small given the relatively low Mach number of the flow). 

However, the main difference between this simulation 
and Model 2a shown in Fig. ^ is the complex structure 
around the galaxy. There is also a dense region at the galaxy 
centre, with ram pressure stripping primarily occurring from 
the outer layers of the galaxy. Complex Kelvin-Helmholtz 
type instabilities develop at the boundary where stripping 
is occurring. The stripping events result in long-dense fila- 
mentary structures trailing the galaxy. A consequence of this 
stripping is the formation of a dense, potentially observable, 
tafl (Section^. 

For Model lb, within the outer bow-shock is a sec- 
ond shock, resulting from the inflow from material from 
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Density Cross — sections Tenaperature Cross — sections 




12 12 

Distance (10^^ cm) Distance (10^^ cm) 

Figure 3. Left panel: The density profiles along the x-axis for the four simulations shown in Figs. |^ and Fig. ^ Shown are results for 
Model lb (solid line), Model 2a (dotted line), Model 2c (dashed line) and Model 3d (dot-dashed line). Note, that the cluster ambient 
density is function of cluster temperature, with the hottest cluster having the highest ambient density. There is a clear difference between 
these simulations, with ES mode simulations (Models 2a and 3d) showing little density enhancements, while MR mode simulations 
(Models lb and 2c) show pronounced density enhancements in both the galaxy region and the downstream region. In particular, the 
double shock structure, discussed by Balsara et al. (1994), is visible for Model lb (solid line). Right Panel: The temperature profiles 
along the x-axis. The lines as the same as in the left panel. The ES mode simulations show only small temperature perturbations, while 
MR mode simulations show more complex structure. 



the downstream region passing through the galaxy (and be- 
ing mass- loaded by the galactic mass-replenishment) before 
interacting with the incident ICM. The density profile for 
Model lb is shown in Fig. ^ and shows both shocks. The 
formation of this second shock is discussed in Balsara et al. 
(1994). 

In Model lb the ram-pressure stripping is sporadic, 
with large amounts of material being stripped over short 
timescales. These mass-stripping events occurs on a time 
scale of a few times the flow-crossing timescale, Tfiow 

To further illustrate the complex flow dynamics for 
Model lb in Fig. ^ we show the x component of the velocity 
at four different locations in the flow. In the undisturbed 
region ahead of the galaxy and bow-shock the gas veloc- 
ity is essentially constant with time. At the galaxy centre 
there is sufficient build up of material to ensure that any 
velocities are small. There are, however, velocity excursions 
with both positive and negative values. The most interesting 
flow pattern is behind the galaxy, where the flow is complex, 
with periods of flow into the galaxy from the downstream 
region. Apparent from Fig. ^ the timescale for these events 
is ~ 5 X 10^ years, a few times Tfiow This is instability 
behaviour noted by Balsara et al. (1994). 

For the case of Model 2c, the basic density and tempera- 
ture pattern is broadly similar to that for Model lb (Fig. |^, 
in that there is a bow-shock and an enhanced density tail. 
However, Model 2c is a less extreme case than Model lb, 
with the enhancement in density at the centre of the galaxy 
being only a factor of a few rather than several orders of 
magnitude. The tail is also narrower and shows less com- 
plex structure. Model 2c seems to be close to the transi- 
tion between the ES and MR mode. At earlier times in the 
simulation it more closely resembles the ES case, but as 
the simulation progressed the velocity at the centre of the 
galaxy dropped with progressive mass-loading. Mass accu- 



mulates at the centre, making it more difficult to strip the 
galaxy (and so on). This simulation illustrates that the dis- 
tinction between the two modes of behaviour is not com- 
pletely clear cut, and in particular as the environment of 
the galaxy changes (as it moves closer or further from the 
galaxy centre) then the dynamics will also change. So, for 
example, while all runs for the hot cluster in Table |^ show 
the ES mode, even with the high mass-replenishment rate, 
if the galaxy were at a larger radius from the cluster centre 
then, because the cluster ram-pressure is lower, the galaxy 
would be able to retain much of the replenished ma ss, and 
the dynamics would be very different (see Section 3.3). It will 
also be possible for the ffow pattern to undergo the oppo- 
site transition. Consider a situation like Model lb, with the 
galaxy losing mass by discontinuous stripping events. If the 
incident ICM ram-pressure were to increase (either higher 
density or higher velocity) , as might be expected if a galaxy 
was falling into the centre of a galaxy cluster, at some point 
the ram pressure would be sufficient to completely strip the 
galaxy, and the galaxy would return to the efficient contin- 
uous stripping mode. Calculations by Takeda et al. (1984) 
suggest that this transition can happen very rapidly, with 
the expulsion of a dense blobs of gas from the centre of the 
galaxy. 



3.3 General Considerations 

Although we have characterised the behaviour of the sim- 
ulations presented here into two broad classes, the ES and 
MR modes, in reality the distinction is not always clear cut, 
as illustrated by Model 2c. 

In Fig. ^ we plot the mass-retention fraction ($) ver- 
sus the ram-pressure stripping parameter (^) as defined by 
Gaetz et al. (1987). The best-fit model for the Gaetz et al. 
(1987) simulations is also shown (eqn.^). The results of our 
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Figure 4. The time-variation of the x-component of the gas ve- 
locity with respect to the galaxy at several points in the flow for 
Model lb. The sampling positions are at a; = 10^^ cm (i.e. well in 
front of the galaxy - dot-dashed line), a; = 9 X 10^'^ cm ( i e. a small 
distance in front of the galaxy centre - dotted line), 10^^ cm (at 
the galaxy centre - dashed line), and 1.5 X 10^^ cm (behind the 
galaxy - solid line). Behind the galaxy the x-component of the 
velocity changes sign several times in a quasi-periodic manner, 
implying back-flow into the galaxy, which results in the build-up 
of material at the centre of the galaxy. This is the 'pumping mode' 
instability discussed in Balsara et al. (1994). 
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Figure 5. The mass retention fraction $ versus the mass- 
stripping parameter ^. The points are the results from the simu- 
lations presented here and the line is the best flt to the results of 
Gaetz et al (1987; eqn.^. 



simulations are similar to those of Gaetz et al. (1987), and 
broadly speaking the transition from the MR mode to the 
ES mode occurs for values of ^ ~ 5 — 10. 

Consideration of Fig. ^ and the simulations presented in 
Balsara et al. (1994) shows that for the MR mode material 
accumulates in the center of the galaxy, and ram-pressure 
stripping occurs from the outer edges of this denser central 




region. So, even while a galaxy may retain mass at its cen- 
tre, beyond a certain radius it will be completely stripped. 
Using the results of Takeda et al. (1984) we can evaluate 
the conditions from ram-pressure stripping at different radii 
within the galaxy. 

Takeda et al. (1984) found that for a spherical galaxy, 
the ISM will be stripped from a galaxy at a projected radius 
d when 



(19) 



where pd is the density of the incident ICM, Vgai is the 
galaxy velocity, a, is the fraction of the stellar mass in- 
jected into the galaxy ISM (i.e. for Mrep = ^ Mq yr~^ , 
a* = 2.6 X 10"'^" s~^), E,(d) is the stellar surface den- 
sity at a projected radius d, a* the stellar velocity disper- 
sion, and C a constant to be empirically determined. For 
Mgai ~ 10^^ Mq a reasonable value is a, = 320 km s~^ 
(c./.Djorgovski & Davis 1987). 

From the stellar distribution used in the hydrodynamic 
model (Section |^ 



E.(d) 



VMTd? 



(20) 



with p*(0) the galaxy central stellar density. In these simu- 
lations a value of C ~ 70 gives a reasonable representation 
of the depth of ram-pressure stripping in the galaxy. For ex- 
ample, for Model lb this corresponds to the galaxy being 
completely stripped beyond a radius of d ~ 7Rc, and for 
Model 2c for d ~ Rc. As discussed by Gaetz et al. (1987), 
this criterion implies that it is energy transfer from the in- 
cident ICM that is important in the stripping process. This 
criterion, although having a slightly different functionality 
than the results from Gaetz et al. (1987), can also be used 
to evaluate the general behaviour of the simulations. If the 
equality in eqn. (|l^) is satisfied for all d then the galaxy 
will be completely stripped. This defines a minimum veloc- 
ity Vstrip, SO that if Vgai > Vstrip then the galaxy will be 
completely stripped. On the other hand, if the equality is 
only satisfied for d > Rh then the galaxy will essentially 
retain all of the replenished mass. 

In Fig. ^ we show the conditions for total stripping 
based on the parameters used in this paper. The relation- 
ship between Vgai and kT^ is shown for the cases when 

Vgai — Vmean and Vgai = 1.5Vmean. AlsO shoWU is Vstrip 

as a function of kTd. Essentially, Vstrip oc T~i'^'^. The case 
when the galaxy is at a radial distance of r = Rcore from 
the cluster centre and Mrep ~ 1 Mq yr~^ is shown (dashed 
line), and suggests that for a cluster with a temperature 
above 3 keV the galaxy will be completely stripped (for 
Vgai ~ Vmean). FoT sl larger mass-replenishment rate of 
Mrep = 5 Mq yT~^ the required temperature for complete 
stripping rises to ~ 5 keV (dotted line). In the simula- 
tions presented here, for the kTci = 4 keV cluster and for 



Mr 



IM, 



yr 



the galaxy is completely ram-pressure 



stripped, while for Mr, 



5 Mq yr the galaxy is stripped 



for the larger value of Vgai but not for the smaller value. 
Moving the galaxy out to larger distances from the clus- 
ter centre increases the stripping temperature (for large r 
Vstrip OC r^^'^), and means that even galaxies in hot clusters 
will not necessarily be completely stripped if they are on the 
periphery of the cluster (dot-dashed line) . It is worth noting 
that for the Fornax and Virgo clusters (with kT ~ 1.3 keV 
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Figure 6. The conditions for total ram-pressure stripping from 
galaxies as a function of the cluster temperature T^;. Shown are 
the galaxy velocities Vg^i as a function of T^j (Section 2.3.2) 
for the two assumed conditions Vgai = Vmean and Vgai = 
l.5Vmean (solid lines). The value of Vgtrip for the case of Mrep = 
1 Mq yr"-"^ is shown (dashed line), as well as the cases when 
Mrep = 5 Mq yr~^ (dotted line). Also shown is the case of 
Mrep = 1 Mq yr"""^ , but with the galaxy at r = bRcore from 
the cluster center instead of r = Rcore (dot-dashed line). For 
complete ram-pressure stripping Vgai > Vstrip- (see Section 3.3 
for more details. 



and 1.9 keV respectively) that for Mrep = ^ Mq yr~^ then 
Vstrip ~ 1680 km and 1350 km respectively. 

The main point from Fig. ^ is that the formation of 
dense tails behind galaxies in clusters is most likely to occur 
in cool clusters, for galaxies with high mass-replenishment 
rates, or alternatively in the outer regions of richer clusters. 

We also note the following general trends from the simu- 
lations. Gas directly behind the bow-shock is always hotter 
than ambient cluster. On the other hand, gas from ram- 
pressure stripping will be cooler than the ambient cluster 
gas {Trep ^ Tci). Depending on the balance between the 
two sources of material the tail will be either hotter than 
(or comparable to) the cluster or cooler than the cluster if 
substantial stripping is occurring. This process is illustrated 
in the values of Ttaii in Table ^, where for MR mode simu- 
lations Ttaii < Tel , while for the ES mode Ttaii ~ Tei . 

One initially counter-intuitive result that follows from 
this is that those galaxies which retain a larger fraction of gas 
(and consequently suffer less from ram-pressure stripping) 
have denser tails, which in turn translates into more visible 
tails. Correspondingly, galaxies that lose the most material 
into the downstream region have less dense tails. The reason 
is straightforward - the density (and visibility) of the down- 
stream region is a sensitive function of the gas-dynamics 
behind the galaxy. Those galaxies that have higher values of 
Mrep have more complex flow patterns behind the galaxy, 
with back- flow occurring and generally lower velocities in the 
tail. This should be contrasted with low mass-replenishment 
cases where the flow is simple - there is some deceleration 
associated with the galaxy, but the a;-component of the ve- 
locity is always positive (c./. Fig. ^). This means that mate- 



rial stripped from the galaxy is moved swiftly away from the 
galaxy and the density never builds up to significant levels. 

As might be expected there is a correlation between the 
opening angle of the bow-shock cone and the Mach num- 
ber. For these simulations we find that the opening angle 
for Mach numbers of ~ 1.3 is around 50°, decreasing to 
around 30° for Mach numbers around 2. While in theory this 
would allow a method for determining the Mach number of 
the galaxy through the ICM (and hence the actual 3-D ve- 
locity of the galaxy) there are some complications. If there 
is sufficient mass-loss from the galaxy then this will alter 
the opening angle. In cases where the mass-replenishment is 
substantial, the cone opening angle is broader, which would 
result in an underestimate of the galaxy Mach number. As 
discussed later, the bow-shock will usually have a faint sig- 
nature and may be difficult to observe. 



4 OBSERVABLE CONSEQUENCES 

One of the main purposes of this paper is to investigate the 
expected observable properties of the theoretical simulations 
of the galaxy /ICM interaction, with an emphasis on the X- 
ray waveband. As most of the current data have been taken 
with the ROSAT satellite (the 0.1 - 2.5 keV waveband) this 
is the spectral region we concentrate on, though future work 
will look at both the XMM and Chandra wavebands. 

The hydrodynamic simulations (Figs. ^ and ^ provide 
values of p and T at all locations, and then it is straightfor- 
ward to generate synthetic X-ray surface brightness maps for 
the simulations. The basic synthetic X-ray surface brightness 
maps for the four simulations in Figs. |l] and ^ are shown in 
Fig. 1^ These images are for the 0.1 — 2.5 keV waveband, and 
the background cluster emission (appropriate for t he w ave- 
band and for the cluster surface brightness - Section 2.2) has 
been added. The cluster emission is often greater than the 
emission from the wake or bow-shock region and this has the 
effect of making the contrast between the tail or wake and 
background emission small. The effects of Galactic absorp- 
tion are not included in Fig. ^ as they are likely to be of little 
consequence. It is necessary to carefully look at the scales in 
Fig. 1^ to determine the visibility of the various structures. 
The images can either have a large (i.e. Model lb) or small 
(the ES mode simulations) dynamic range, and the scaling 
has been adjusted appropriately. The important features to 
notice in Fig. ^ are as follows: 

(i) For the ES mode simulations the most prominent 
structures are the bow-shocks, with a small deficit of emis- 
sion preceding the bow shock. This is due to the lower den- 
sity region caused by the gravitational acceleration of mate- 
rial by the galaxy. 

(ii) For the hotter cluster simulations (see Model 3d), the 
cluster background is very bright, and the contrast between 
bow-shock and background is small (~ 1 per cent). This 
is more apparent in Fig. ^ where surface brightness slices 
are shown for all 4 simulations. This means that interaction 
effects will be more likely to be seen in cooler, less luminous 
clusters. 

(iii) For the MR mode (Models lb and 2c) there are three 
main features visible in the synthetic images. First, the bow- 
shock, second elongated emission from the tail region and 
third emission from the galaxy itself. 
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Figure 7. Synthetic ROSAT surface brightness profiles for the four simulations shown in Figs, and^. Top left panel: Model 2a; top 
right panel: Model 3d (both ES mode). Bottom left panel: Model lb; bottom right panel: Model 2c (both MR mode). The synthetic 
surface brightness images have been calculated in the 0.1 — 2.5 keV waveband, with units of erg cm~^s~^arcmin~^ . Note that emission 
from the cluster has been included. Careful attention should be paid to the scaling of these images, as some of them have a much larger 
dynamic range than others. 



(iv) For Model lb, the galaxy is very luminous, and in or- 
der to enhance the visibility of the tail and bow-shock region 
we do not show the full dynamic range (c./. Fig. |^. There is 
a large accumulation of material at the galaxy centre. How- 
ever, this emission, essentially analogous to a cluster cooling 
flow, will likely be strongly absorbed in a real galaxy. 

(v) For Model 2c, both the bow-shock and tail are visible, 
but the tail is less pronounced that for Model lb. 

(vi) For both Model lb and Model 2c the emission from 
the tail does stand out substantially from the cluster back- 
ground. However, the contrast for the bow-shock is less. This 
seems to be a general characteristic - the bow-shock is al- 
ways less visible than any ram-pressure stripped tail. 

(vii) For Model lb the synthetic image shows considerable 
substructure on a scale of around 5 — 10". Such structures 
should be clearly visible with Chandra and XMM. 

To further illustrate these points, in Fig. ^ surface 
brightness slices are shown. These slices have been gener- 
ated by taking a slice of width 10^^ cm along the a;-axis. For 
the ES cases (Models 2a and 3d as examples) the increase 
in surface brightness is relatively small (a few per cent at 
most) and the surface brightness profile is relatively simple. 
For the MR mode case there are much greater increases in 
emission. Taking Model lb as an example, the region associ- 
ated with the galaxy is extremely bright, but the tail region 



also shows a substantial increase, with an excess brightness 
over the background of ~ 10 — 30 per cent. Even the case of 
Model 2c, which shows less in the way of enhanced emission 
at the position of the galaxy, has a tail for which the surface 
brightness is more than the ~ 10 per cent greater than the 
underlying cluster. 

For both the MR mode simulations in Fig. ^ there is 
substantial substructure in the surface brightness in the tail 
region, with regions of enhanced emission associated with 
blobs of gas stripped from the galaxy. The bow-shock region 
in the MR mode cases extends further in front of the galaxy 
than for the ES mode cases. 

For our simulations we have generated hardness ratio 
maps, with the hardness ratio defined as {H — S)/{H + S), 
where H is the hard-band (0.5 — 2.5 keV) image and S the 
soft-band image (0.1 — 0.5 keV). The results in Table |^ show 
that for the MR mode the tail region is cooler than the am- 
bient cluster, and consequently spatial hardness variations 
might be expected. However, because the temperature dif- 
ferences are usually small and the emission from the tail 
region is always small compared to the background cluster 
emission, this means that differences in the hardness ratio 
are always small. 
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Figure 8. Surface brightness slices for tlie four images shown in Fig. ^ The units are erg cm~^s~^arcmin~^ . These slices in surface 
brightness are taken along the x-axis through the galaxy for a slice with a width of 10^-^ cm. The level of background cluster emission 
varies substantially from the cool to the hot cluster. Also note that the scale for Model lb has been truncated to enhance the contrast in 
the tail region. The peak surface brightness for this model is ~ 10~^^ erg cm~-^s~^arcmin~^ . Note the substructure in the tails of the 
MR mode simulations. 



5 DISCUSSION 

In this paper we have presented some new simulations of the 
interaction between galaxies and the ICM in clusters. For the 
first time we have attempted to quantify the likely observa- 
tional characteristics of structures resulting from such in- 
teractions, with some promising results. Galaxies which are 
capable of retaining some mass and have substantial tails are 
likely to have features that are observable. However, those 
that are completely ram-pressure stripped are unlikely to 
have observable features. The bow-shock is a more promi- 
nent feature in MR mode galaxies than in ES mode galaxies 
(Figs. I and I). 

The likely criterion for galaxies to show observable in- 
teraction effects in clusters are as follows: 

(i) Galaxies with younger stellar populations (and hence 
more substantial mass-loss) are more likely to show interac- 
tion effects at X-ray energies. 

(ii) Galaxies in cooler clusters or group will be more likely 
to have observable stripped tails. Galaxies in rich clusters 
will be more likely to be completely stripped. 

(iii) Galaxies in lower surface brightness clusters - in 
richer, brighter clusters the contrast between tail and back- 
ground emission is smaller. 



(iv) Galaxies at the periphery of rich clusters. Massive 
gas-rich galaxies infalling into rich clusters should show tails. 

This means that bluer galaxies in nearby cooler clusters 
are the most likely place to look for such interaction effects. 
For rich clusters, such as Coma, most galaxies are likely to 
be completely stripped, except those at the periphery of the 
cluster, which may be in an environment where they can 
retain some mass. These simulations have also shown that 
while tails may be a feature of galaxies in clusters of all 
levels of richness the effects are more likely to be seen in the 
cooler clusters. 

The potential rewards of being able to exploit the 
galaxy/ICM interactions are great (Merrifield 1998), how- 
ever these simulations show the dynamics are complicated 
and the level of emission from the tail (for instance) will de- 
pend on parameters that may be difficult to determine (such 
as Mrep). Consequently, it may take a substantial effort at 
modelling individual galaxies to fully determine the true ve- 
locity vector for that galaxy. The other, simple approach, to 
merely use the galaxies direction of motion on the plane of 
the sky (as revealed by a linear trailing structure), may also 
be able to yield interesting results. 

The assumptions made here have been chosen to be 
simple and to be the most favourable. For example, we have 
assumed the galaxies to be spherical. In reality, the galax- 
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ies will be elliptical and the gas flow out of the galaxy will 
depend on the orientation of the galaxy with respect to its 
velocity vector through the ICM. For instance, it might be 
anticipated that mass might be lost more easily from the 
galaxy through the minor axis of the galaxy. We have also 
assumed that we are viewing the galaxy in the plane of the 
sky. If the galaxy is moving at a different angle then projec- 
tion effects will play a role, foreshortening the tail and po- 
tentially making identification of galaxy/ICM features dif- 
ficult. The effect of both of these assumptions will have to 
be explored with more sophisticated simulations. While cer- 
tain orientations of an elliptical galaxy can be dealt with 
by 2-D simulations, the more general case of random galaxy 
orientation will require 3-D simulations 

For several of the suspected examples of galaxy/ICM 
interaction (Section |l|) all that can be discerned about the 
dynamic structures associated with the galaxy motion is an 
extension of the X-ray contours away from the cluster centre 
(usually interpreted as due to infall of the galaxy towards 
the cluster centre). Examples of this are NGC4839 in the 
Coma cluster (Dow & White 1995) and NGC 1404 (Jones 
et al 1997). It is notable that NGC 4839 is located towards 
the edge of the Coma cluster, where we might expect inter- 
actions effects to be more visible. 

In the case of M86, in addition to an emission peak 
at the galaxy centre, there is also a second peak associated 
with the tail. These simulations and those of Balsara et al. 
(1994) suggest that stripping will often occur in discrete 
events, with blobs of material being removed. These blobs 
will result in enhancements in the tail region, and this is 
what may be occurring in M86. 

NGC 4472 shows very complex emission, with some ev- 
idence of a bow-shock as well as a tail. This galaxy will 
clearly require more detailed work on it. Irwin & Sarazin 
(1996) suggest that some of the features may be related to 
the elliptical nature of the galaxy and its orientation with 
respect to its motion. This idea will clearly have to be tested 
with further simulations allowing for non-spherical galaxies. 

In summary, we have new presented calculations of 
the galaxy/ICM interaction in clusters, with particular em- 
phasis on investigating the likely observable consequences 
of such interactions. From these simulations it is clear 
that many galaxies should show observable effects of the 
galaxy/ICM interaction, the most prominent feature being 
linear tails, resulting from ram-pressure stripping. We have 
performed simulations for a wide range of cluster param- 
eters, with the goal of trying to define the best parame- 
ter space to investigate. While massive galaxies, with young 
stellar populations, in cool nearby clusters such as Virgo or 
Fornax offer the best possibilities, the periphery of richer 
clusters also offers interesting possibilities. An important 
result is that it is not galaxies that are completely ram- 
pressure stripped that have the brightest tails, but rather 
those which are able to retain a substantial fraction (but not 
all) of the replenished material. In this situation, stripping 
occurs in discrete events and results in high density tails. 
Bow-shocks should also be visible in some galaxies, particu- 
larly in upcoming XMM and Chandra observations of nearby 
candidates. There are good prospects that, through a com- 
bination of modelling and observations with the new X-ray 
observatories, galaxy/ICM interactions offer a new window 
to investigating cluster dynamics. 
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